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ABSTRACT 

Recent observational studies of w Centauri by Hubble Space Telescope have discovered a double main 
sequence in the color magnitude diagrams (CMDs) of its stellar populations. The stellar population 
with the blue main sequence (bMS) is observationally suggested to have a helium abundance much 
larger, by AY ^ 0.12, than that of the red main sequence (rMS). By using somewhat idealized models 
in which stars of the bMS are formed from gas ejected from those of the rMS, we quantitatively 
investigate whether the helium overabundance of the bMS can result from self-enrichment from massive 
AGB stars, from mass loss of very massive young stars, or from type H supernovae within w Gen. We 
show that as long as the helium enrichment is due to ejecta from the rMS formed earlier than the 
bMS, none of the above three enrichment scenarios can explain the observed properties of the bMS 
self-consistently for reasonable IMFs. The common, serious problem in all cases is that the observed 
number fraction of the bMS can not be explained without assuming unusually top-heavy IMFs. This 
failure of the self-enrichment scenarios implies that most of the helium-enriched gas necessary for the 
formation of the bMS originated from other external sources. We thus suggest a new scenario that 
most of the second generation of stars (i.e., the bMS) in u) Gen could be formed from gas ejected from 
field stellar populations that surrounded uj Cen when it was a nucleus of an ancient dwarf galaxy. 

Subject headings: globular clusters: general - globular clusters:individual (a; Gentauri) 


1. INTRODUCTION 

One of the most remarkable results of recent observa¬ 
tional studies of oj Gen is that it shows a double main se¬ 
quence (DMS) in the color magnitude diagrams (GMDs) 
of its stellar populations (e.g., Anderson 1997; Bedin et 
al. 2004). Bedin et al. (2004) proposed four possi¬ 
ble scenarios for the origin of the DMS: (1) Theoreti¬ 
cal isochrone models or data calibration are in error, ( 2 ) 
Stars on the bluer main sequence (bMS) of the DMS are 
a super metal-poor ([Fe/H]^ —2.0) population, (3) The 
bMS represents a very helium-rich (Y > 0.3) population, 
and (4) The bMS represents a background stellar popu¬ 
lation about 1 — 2 kpc behind uj Gen. A number of recent 
investigations have suggested that the above scenario (3) 
(hereafter referred to as the “helium pollution scenario” 
and abbreviated to “HEPS” for convenience) is the most 
promising among the four (e.g., Lee et al. 2005; Norris 
2004; Piotto et al. 2005). 

In the HEPS, there were two epochs of major star for¬ 
mation in the history of oj Gen. The first generation 
of stars was formed from proto-globular cluster (GG) 
cloud(s) to become stars on the red main sequence (rMS) 
whereas the second generation of stars was formed from 
gas ejected from the rMS. One of the key questions re¬ 
lated to the HEPS is whether and how the observation- 
ally suggested large helium enrichment (AY 0.12) for 
the bMS can be obtained from the rMS with normal Y 
(~ 0.24) in the star formation history of w Gen (e.g.. 
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Norris 2004). There are three candidates for the helium 
overabundance (Norris 2004; D’Antona et al. 2005; Pi¬ 
otto et al. 2005): (1) AGB stars with the masses larger 
than 6 M 0 , ( 2 ) stellar winds associated with massive stars 
during their early evolutionary phases, and (3) type H 
supernovae (SNe H). It is however unclear which of the 
above three is the most reasonable and realistic in the 
HEPS given a lack of extensive theoretical studies on the 
three candidates. 

The purpose of this Letter is to investigate the three 
possibilities in a quantitative manner and thereby dis¬ 
cuss which is the most promising as the cause of the 
bMS (and the DMS) in the context of the HEPS. In this 
investigation, (1) the possible Y value observationally 
suggested for to Gen (Norris 2004; Piotto et al. 2005) 
and (2) the observed number fraction of the bMS (Be¬ 
din et al. 2004) are used to constrain the best possi¬ 
ble initial mass function (IMF) of forming stars in the 
HEPS. We do not intend to discuss extensively the ob¬ 
served abundance pattern of the bMS and the rMS in this 
paper, because chemical yield tables for AGB stars with 
helium-rich ejecta (Y > 0.35) are not currently avail¬ 
able. We will discuss this point in our forthcoming pa¬ 
pers (Bekki & Norris 2005; BN). Previous theoretical 
studies demonstrated that if GGs lose more than 50% of 
their initial masses, they will disintegrate (e.g., Geyer & 
Burkert 2001). We also use this result as a constraint for 
globular cluster IMF in the HEPS. We demonstrate that 
none of the above three candidates can explain the above 
constraint (2) without the modeled uj Cen disintegrating. 
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2. THE POSSIBLE THREE POLLUTERS 

Firstly we summarize the possible three “polluters” 
which can significantly enrich u) Cen with helium if the 
second generation of stars (i.e., the bMS) was formed 
from their ejecta. In order for the observationally sug¬ 
gested large value oiY (~ 0.35) to be explained, we need 
to select stellar objects whose ejecta have Y significantly 
higher than 0.24. Figure 1 shows the helium mass frac¬ 
tion (Y) in the ejecta of three candidates; (1) AGB stars 
(van den Hoek & Groenewegen 1997, hereafter referred 
to as VG97); Ventura & D’Antona 2005, as VA05), (2) 
massive stars (Schaller et al. 1992, S92), and (3) SNe 
II (Woosley & Weaver 1995, WW95). Based on the val¬ 
ues listed in these papers, Y is derived for each stellar 
mass. The theoretical predictions of Y for AGB stars is 
at most 0.28 for VG97 and 0.32 for VA05. D’Antona et 
al. (2005) however suggested that more massive AGB 
stars with the masses of 6 — 7M0 can have Y = 0.40 (see 
also Lattanzio et al. 2004). We therefore consider that 
massive AGB stars can be the polluters in oj Gen but 
that lower mass AGB objects (< 5 M 0 ) with smaller Y 
can not. 

The observed difference in metallicity ([Fe/H]) between 
the rMS and the bMS (Piotto et al. 2005) suggests that 
the above scenario (1) alone can not simply explain the 
metallicity difference, because metallicity of AGB ejecta 
do not differ from those of their host stars. Therefore, 
full chemical evolution models including the contribution 
of the three polluters (and mixing between the ejecta 
of the polluters and fresh gas) need to be considered if 
we try to explain fully self-consistently both the helium 
and metallicity one between the rMS and the bMS. We 
however discuss advantages and disadvantages of each 
of the three polluters in the helium enrichment process 
separately. This will better permit us to disentangle the 
contribution of the different polluters in producing He- 
enhanced stars and thereby to analyze advantages and 
disadvantages of each scenario in explaining the helium 
abundance of the bMS. The entire chemical evolution 
history for different elements (e.g., ^^G) in uj Cen will be 
discussed in our future papers (BN). 

2.1. AGB stars 

The first question is whether the number fraction of 
the second generation of stars (i.e., the bMS fraction) can 
be as large as 0.25 (Bedin et al. 2004) for a reasonable 
IMF without OJ Cen being disintegrated following mass 
loss from massive stars and SNe II. In order to estimate 
the mass fraction (/agb) of AGB progenitor stars with 
masses ranging from 6 M 0 to 7 M 0 in a GC with the total 
mass of Mci, we assume an IMF in number defined as 
tp(mi) = where mi is the initial mass of each 

individual star and the slope s = 2.35 corresponds to the 
Salpeter IMF. The normalization factor A is a function 
of Mci, mi (lower mass cut-off), and (upper one). 
A is expressed as A = ' "''’here m\ and mu 

are set to be O.IM 0 and I 2 OM 0 , respectively. Although 
the number fraction of low-mass stars on the bMS can 
depend on the forms of IMFs, we adopt the above IMF 
and show the results. This is mainly because we can show 
more clearly the roles of IMF slopes in controlling the 
mass fraction of AGB stars for the adopted IMF in this 
preliminary study. The models with the IMF proposed 


by Kroupa et al. (1993) show a rather small fraction 
(~ 0.1 %) of AGB stars in comparsion with the models 
with the above IMF. 

In this AGB pollution scenario, we assume that (1) the 
gas ejected from AGB stars in the first generation (i.e., 
the rMS) is all used (i.e., a star formation efficiency of 
100 %) for the formation of the second generation (i.e., 
the bMS) and (2) the total amount of gas ejected from 
a AGB star is the same as the stellar mass. This second 
assumption is reasonable because more than 80 % of the 
stellar mass is ejected as AGB winds for AGB stars with 
mi > 6 M 0 (e.g., VG97). We also assume two different 
slopes of IMF (s) for the first (si) and the second (S 2 ) 
generations of stars in order to investigate the maximum 
value of / 2 nd/(/ist + / 2 nd), where /ist and / 2 nd are the 
number of stars (per unit mass) with masses smaller than 
O. 88 M 0 (corresponding to stars older than ~ 12 Gyr) for 
the first generation and for the second one, respectively. 

Figure 2 shows that the smaller the IMF slope si is 
(i.e., the IMF becomes more “top-heavy”), the larger the 
number fraction of the second generation (/ 2 nd/(/ist + 
/ 2 nd)) becomes for a fixed S 2 of 2.35. The ratio of the 
final cluster mass to its initial mass (/rem) however also 
becomes smaller for smaller si. As a result of this, /rem 
of the models with si < 1.75 become smaller than the 
threshold value (0.5) below which star clusters can be 
disintegrated as the result of mass loss (e.g., Geyer & 
Burkert 2001): No models can be located within the re¬ 
gions with 0.25 < / 2 nd/(/ist + / 2 nd) < 0.35 and with 
/rem > 0.5 where oj Gen can survive possible disintegra¬ 
tion and have the bMS fraction consistent with obser¬ 
vations. The number fraction of the second generation 
can not be larger than 0.16, even for §2 = 5.35 (i.e., ex¬ 
tremely bottom-heavy) and si = 1.95 for which oo Cen 
can manage to survive disintegration due to mass loss. 
Thus the results in Figure 2 suggest that the observed 
fraction of the bMS can not be explained simply by the 
AGB pollution scenario for reasonable IMFs. 

If ruu is as low as 8 M 0 , oj Cen can not disintegrate 
owing to mass loss from SNe II for any values of si and 
S 2 - It is found that / 2 nd/(/ist + / 2 nd) is about 0.29 for 
the model with TOu = 8 M 0 , si = 1.35, and S 2 = 2.35. 
This result suggests that if oo Cen was formed with a very 
unusual IMF (i.e., with no or little SNe II), the observed 
bMS fraction can be reproduced without disintegration 
of OJ Cen. The above model with no or little stellar rem¬ 
nants (black holes and neutron stars) of massive stars, 
however, would not be consistent with the presence of 
neutron stars in oj Cen (e.g., Haggard et al. 2004). 

2 .2. Massive stars 

The helium mass fraction in the stellar winds of mas¬ 
sive stars with mi > 85M0 can become as high as 0.4 
with a maximum of 0.49 for mi = I 2 OM 0 (Schaller et 
al. 1992). Although there is no problem with Y in 
this scenario, the mass fraction of massive stars with 
mi > 85M0 can not be as large as ~ 0.25 for reason¬ 
able IMFs (1.95 < s < 2.35) for which oj Cen does not 
disintegrate. Figure 3 shows that the mass fraction of the 
ejecta of the massive stars is at most 0.2 for a plausible 
range of IMF slopes. Given the fact that the observed 
star formation efficiency in the Galaxy can be at most 
~ 0.4 (Wilking & Lada 1983 for the p Oph cloud), this 
result suggests that the observed fraction of the bMS 
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(~ 0.25) is unlikely to be explained by this scenario. 

A more serious problem for this stellar wind scenario 
is that there would be no or little difference in the mean 
metallicity between the first generation of stars and the 
ejecta of these objects (Schaller et al. 1992): Piotto et al. 
(2005) demonstrated that there is a metallicity difference 
by A[M/H]~ 0.3 between the rMS and the (more metal- 
rich) bMS in uj Cen. Another possible problem of this 
stellar wind scenario is that the ejecta should be very 
quickly converted into stars (within less than ~ 10^ yr) 
between the stellar wind phases of very massive stars and 
the onset of SNe II that can blow away the ejecta from 
Lu Cen. 

2.3. Type II supernovae 

The helium mass fraction (T) in the ejecta of SNe II 
stars with 12 < toi/M© < 40 for metal-poor stars with 
Z = 0.01^0 ranges from 0.37 to 0.45. (WW95). Al¬ 
though this scenario has no problems with Y, the abun¬ 
dance by mass for heavier elements in the ejecta is too 
high to be consistent with observations. and ^°Ca 
in WW95 are as high as ~ 0.01 and 0.001, respectively, 
which are higher than the solar values (~ 3.0 x 10“^ 
and ~ 6.0 x 10“^, respectively; WW95) and thus much 
higher than those of u Cen. This inconsistency with the 
observed metallicities suggests that the bMS can not be 
formed directly from the ejecta of SNe II of the rMS. If 
the ejecta of SNe II can be mixed with the primordial 
gas used for the formation of the rMS, the abundance of 
’^^C and '^°Ca of the bMS formed from the mixed gas can 
be significantly reduced. However, the helium abundance 
would also be significantly reduced in the process of mix¬ 
ing. Therefore it can be concluded that this scenario is 
highly unlikely to explain the observed abundance of the 
bMS. Furthermore this scenario has difficulty in explain¬ 
ing the observed fraction of the bMS without assuming 
top-heavy IMFs, just as the other two scenarios do. 

3. DISCUSSION: ALTERNATIVE SCENARIOS 
3.1. External pollution in the nucleus of a dwarf galaxy 

The present study shows that the proposed three sce¬ 
narios all fail to explain self-consistently the observed 
properties of the bMS in uj Cen. The AGB scenario ap¬ 
pears to be the most promising among the three, though 
it remains unclear whether this scenario can explain the 
observed abundances other than helium in to Cen. The 
most serious problem of the AGB scenario is that the 
total mass of the ejecta of more massive AGB stars from 
the rMS is too small to be consistent with the observed 
fraction of the bMS of uj Cen. If however we relax the 
adopted assumption that all stars of the bMS originate 
from AGB ejecta of the rMS initially within oj Cen, the 
problem of the AGB scenario can be significantly allevi¬ 
ated. The problem of the bMS fraction can be solved if 
the original total mass of the rMS was larger (by a fac¬ 
tor of ~ 10) than the present mass (i.e., if uj Gen was a 
super-giant GG) and if most of the stars of the rMS were 
removed from to Cen after their ejecta were used for the 
formation of the bMS. 

We here discuss a scenario in which uj Cen was the nu¬ 
cleus of a nucleated dwarf galaxy (dE,N) where the AGB 
ejecta of the central field stellar populations surrounding 
the (proto-) uj Cen were consumed for star formation of 


the bMS. In this scenario, the ancient host dE,N of to 
Cen has been completely destroyed by the Galactic tidal 
field so that only its nucleus (i.e., uj Cen) is now observed 
(e.g., Bekki & Freeman 2003). Not only the ejecta of the 
rMS but also those from stellar populations surrounding 
the nucleus can be used for the formation of the bMS: 
the above bMS fraction problem should not be so serious 
in this scenario. The bMS can be more metal-rich than 
the rMS if the mean abundance of the field stellar popu¬ 
lations and gas in the dwarf is slightly higher than that 
of the rMS. 

The key question in this scenario is how much gas 
ejected from the central field stars can be converted into 
the bMS without being removed by energetic outflow of 
AGB stars in the central regions of uj Cen’s host dE,N. 
Recent numerical simulations on transformation from a 
dE,N into uj Cen have suggested that w Cen’s host dE,N 
has Mb ~ —14 mag, corresponding to the total stellar 
mass of 1.25X 10® Mq (e.g., Bekki & Freeman 1993). If 
the total mass of the bMS is about 30% of present-day uj 
Cen’s mass (= 5.0 x 10® Mq; Meylan et al. 1995), it is 
necessary for only 1.2% of the stellar mass of the dE,N to 
have been converted into the bMS. Since the mass frac¬ 
tion of more massive AGB stars in the dE,N (> 0.02) 
for a reasonable IMF can be significantly larger than the 
above value (~ 0.01), the observed bMS fraction is not a 
problem in this new scenario. 

The central escape velocity (I4sc) of this host embed¬ 
ded in a massive dark matter halo with the total mass of 
10® — IO^Mq is 50 — 90 km s“^ depending on the inner 
profile of dark matter halo (Bekki 2005). The derived 
I4sc is significantly larger than the observed maximum 
wind velocity of ~ 40 km s“^ for AGB stars (e.g., Loup 
et al. 1993). Therefore, the AGB ejecta of the dE,N 
is highly likely to be trapped within the central region 
of the dE,N and consequently converted into new stars 
(i.e., the bMS). Although extensive numerical simula¬ 
tions on the central gas dynamics of a dE,N are necessary 
to confirm that gas ejected from field stars of the dE,N 
can be transferred into the nuclear region and converted 
into new stars there, we propose that the bMS can origi¬ 
nate from gas ejected from stars within the ancient dE,N 
whose nucleus was uj Cen. 

3.2. Helium sedimentation 

The local helium abundance of the remaining gas of 
forming uj Cen might have been significantly higher (lead¬ 
ing to a bMS with higher helium abundance), if helium 
sedimentation due to gravitational diffusion (which was 
originally proposed by Fabian & Pringle (1977) for a 
mechanism responsible for gaseous abundance gradients 
in clusters of galaxies) occurred in uj Cen. Chuzhoyfe 
Loeb (2004) showed that gravitational diffusion in the 
interstellar medium of giant elliptical galaxies can re¬ 
sult in the increase of helium abundance by a factor of 
1 -I- 0.2(T/10 ^)^’^/Fb, where T and Fb are the gas tem¬ 
perature and the suppression factor (likely to be > 5) of 
the gravitational diffusion by the magnetic field. These 
results suggest that a significant increase of helium abun¬ 
dance is highly unlikely within GCs, which have the cen¬ 
tral velocity dispersions of order of 10 km s“^, corre¬ 
sponding to T ~ 10000 K. 

4. CONCLUSION 
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We have shown that the observed bMS fraction of w 
Cen can not be simply explained by any HEPS in which 
the bMS was formed from ejecta of the rMS. We accord¬ 
ingly have suggested an “external pollution” scenario in 
which the bMS was formed from gas that was initially 
within the central region of ui Cen’s host galaxy. The 
question yet to be answered in this scenario is how star 
formation could proceed within the rMS when gas was 
transferred to the central region. It would be possible 
that the bMS was formed outside yet close to the rMS 
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Fig. 1.— The mass fraction of helium ^He (Y) in the ejecta of stellar winds of massive stars (thick solid), SNe II (thick dotted), and 
AGB stars (thin solid and dotted) for different masses (Mg). For comparison, the observationally suggested value of Y for the bMS of to 
Cen is shown by a dashed line. Y values were calculated by using the tables by Schaller et al. (1992) for massive stars, Woosley Sz Weaver 
(1995) for SNe II, and van den Hoek &; Groenewegen (1997) (thin solid line) and Ventura & D’Antona (2005) (thin dotted line) for AGB 
stars. 
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Fig. 2. — The locations of uj Cen models with different IMFs in the /rem-/ 2 nd/(/ist +/2nd) plane. Here /rem, /ist, and / 2 nd represent the 
mass fraction of remnant stars (i.e., low-mass, long-lived stars) after mass loss due to supernovae, the number of stars of the first generation, 
and that of the second. The vertical dotted line represents the boundary (/rem = 0.5) beyond which lu Cen can survive disintegration due 
to mass loss. The horizontal dotted line represents the lower limit of the observed bMS fraction of lo Cen in Bedin et al. (2004). The upper 
right region surrounded by these two lines is where any uj Cen model should lie to explain observations (i.e., self-bounded and larger bMS 
fraction). The models with si = 1.35 (left), 1.55 (second from left), 1.75 (third from left), 1.95 (third from right), 2.15 (second from right), 
and 2.35 (right) for a fixed S 2 of 2.35 (the Salpeter IMF) are shown by a thick solid line. The models with S 2 = 5.35 (top), 4.35 (second 
from top), 3.35 (third from top), 2.85 (third from bottom), 2.35 (second from bottom), and 1.95 (bottom) for a fixed si of 1.95 are shown 
by a thin solid line. si ( 52 ) is the slope of the IMF in the form of 'i/;(mi) = Ami~^ for the first (second) generation of stars. Note that all 
models have serious difficulties in reproducing both /rem > 0.5 and / 2 nd/(/ist + / 2 nd) ^ 0.25 simultaneously. 
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Fig. 3.— Dependences of the mass fraction of ejecta from massive stars with 85 < mi/M 0 < 120 on the IMF slope (si) of the first 
generation of stars for the models with mi = O.IMq (thick) and mi = 0.8 Mq (thin). 




